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ABSTRACT 


(5-sunspots, with highly complex magnetic structures, are very productive in ener¬ 
getic eruptive events, such as X-class flares and homologous eruptions. We here study 
the formation of such complex magnetic structures by numerical simulations of mag¬ 
netic flux emergence from the convection zone into the corona in an active-region-scale 
domain. In our simulation, two pairs of bipolar sunspots form on the surface, origi¬ 
nating from tw o buoyant segments of a single subsurface twisted flux rope, following 
the approach of lToriumi et al.l (j2014l ). Expansion and rotation of the emerging fields in 
the two bipoles drive the two opposite polarities into each other with apparent rotat¬ 
ing motion, producing a compact (5-sunspot with a sharp polarity inversion line. The 
formation of the (5-sunspot in such a realistic-scale domain produces emerging patterns 
similar to those formed in observations, e.g. the inverted polarity against Hale’s law, 
the curvilinear motion of the spot, strong transverse field with highly sheared magnetic 
and velocity fields at the PIL. Strong current builds up at the PIL, giving rise to recon¬ 
nection, which produces a complex coronal magnetic connectivity with non-potential 
fields in the (5-spot overlaid by more relaxed fields connecting the two polarities at the 
two ends. 


Subject headings: magnetohydrodynamics(MHD) — Sun: sunspots — Sun: photo¬ 
sphere 


1. Introduction 


(5-sunspots, where two spots of opposite polarities share the same penumbrae, are of partic- 
ular interest due to their high productivity of powerful f l ares and coronal mass ejections fCMEs ) 
( Hagvard et al. 1984 : Tanaka 1991 : Schmieder et al. 1994 : Sammis et al. 200(1 : Liu fc Zhang 2001 ). 
Based on their chara cteristics during the formation observed at Big Bear Solar Observatory, 
Zirin &: Liggetti (jl987l l categorize the (5-sunspots into three types: island 6 which emerge at once 
but with different dipoles intertwined, single 5 with large satellite dipoles emerging in the penum¬ 
bra of a large spot, and a third and most popular type formed by collision of two bipolar spot 
groups. Regardless of their types, all the 5 spots, which form from opposite polarities of different 
dipoles and die together, are short-lived complex structures with highly sheared field lines at the 
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PIL. They manifests their complexity on the photosphere as highly structured active regions, with 
non-potential coronal fieMs in the form of sheared and stretched fie ld lines connecting the sunspots 
([Patty Sz Hagvardlll986l : lTanakalll99ll : iLites et alJll995l : lLekalll997l l. storing a large amount of free 
magnetic energy available for the flare activities. In many o ccasions, the footpoints of the d-spots 
exhibit shearing motion at the PIL and rotation of the spots (|Tanakalll99ll : lLeka et al.lll996l ). which 
further stretch and twist the magnetic field lines. The combination of non-potential magnetic field 
and shearing, rotating motion at the PIL produces a magnetic configuration that is prone to mag- 
netic reconnection and eruption, which might be repetitive in the presence of energizing motions 
( Fang et al. 2012a). 


According to the observed features of compact and complex d-spots regions, IZirin Sz Liggett 


(| 19871 ) point out that they may be the result of a subsurface phenomenon which locks the two spots 
closely together and it is therefore important to investigate the sub surface process es tha .t form and 


maint ain these complex spots and build up the free energy. Later iTanakal (| 199 il l and iLeka et al 


(| 19901 ) suggest that 5-sunspots can be formed by the emergence of highly twisted magnetic flux 
tubes from the convection zone, which can explain self-consistently the shearing and rotation mo¬ 
tion as well as the compact pattern of opposite polarities. This concept then motivated series of 
numerical si mulations on helical kink instability of h ighly twisted rising flux tubes emerging from 


the interior ([Linton et al. 


1998 


19991 : 


Fan et al. 


1999|). The kinked evolution of the rising tubes pro¬ 


duces a compact magnetic bipolar region, with polarity order reversed from Hale’s law, with great 
similarities to observed features of 5-sunspots. The braided kinked magnetic stru cture keeps the 
two opposite polarities closely together and drives the rotation around each other. iTanakal ([l99ll ) 
also proposed a model which involved the successive emergences of two twisted magnetic flux ropes 
which may or may not be interconnected under the surface, to explain the evolution of complex 
5-spo ts. Recent study on the quadrupolar active region 11158 with colliding 5-spots ([Sun et al 


20121 ) found that the magnetic fields were highly sheared at the PIL carrying s t rong c urrent, and 
significant free energy was injected into the atmosphere. Later iToriumi et al.l ([20141 1 considered 
two scenarios of the subsurface emerging structures for simulating the formation of 5-spots: two 
flux ropes each with one buoyant section, and a single flux rope with two buoyant sections. They 
found that only the emergence of the two buoyant sections from a single flux rope forms a sharp 
PIL in the domain. Following this work, here we also carry out a simulation of the emergence of a 
twisted flux rope with two buoyantly rising sections, but in a domain of size comparable to active 
regions to study whether compressed 5-spot like structures of realistic active region size scales can 
still form. In the following sections, we describe the numerical setup of our simulation in Section [21 
and investigate the formation of 5-sunpots and the photospheric flows and energy flux in Section [3] 
and m respectively, followed by conclusions in Section [5l 
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2. Simulation Setup 


We solve the ideal MHD equations within the Block Adaptive Tree S olar-wind Roe Up¬ 
wind Scheme (B ATS-R-US) that was developed at the University of Michigan (jPowell et alJll999l : 
Toth et ah 2012): 
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where p, u, E, p, T, B and g are the density, velocity, total energy density, pressure, temperature, 

magnetic held and gravitational acceleration, respectively. The domain of our simulation covers 

158(X)X 138(Y)x79(Z) Mm^, extending 20 Mm below the photosphere and 59 Mm into the c orona 

with an initial backgronnd stratihcation of atmosphere in hydrostatic equilibrium (jFan|[200l|). The 

cell sizes range from 38.6 to 618.68 km in the vertical direction from the interior to coronal height 

and 77.3 to 1237.4 km in horizontal directions. A horizontal hux rope in the x direction is inserted 

2/2 2 / 2 ^ 

at depth Z = -10 Mm, with the magnetic held given by B = BQe~^ x + qrBQe~'' 6, where 
Bq = —12 kG, is the strength of the magnetic helds at the central axis of the rope, 6 is the 
azimnthal direction in the cross section of the hux rope, q = 0.50/a, is the rate of twist per unit 
length, a = 3.0 Mm is the radius of Gaussian decay, and r is the distance from the axis of the hux 
rope. Fig. [T] shows the vertical stratihcation of pressure and temperature in the domain. The green 
line represents the distribution of magnetic energy inside the hux rope in a cross-section cut. The 
plasma beta at the central axis of the hnx rope, /3 = 66. The initial rope is density-depleted at two 
segments centered at x = ±15 Mm as shown in Fig. O with perturbation as: 
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otherwise. 


( 5 ) 


Here po is the initial density value in hydrostatic equilibrium, buoyancy factor p = 0.15, and x± 
= ± 15 Mm, A = 4.6 Mm is the half length of each buoyant segments. And the thermal pressure 
inside the hux rope is modihed to maintain the total pressure balance. 


3. Formation of 5— sunspots 

The density-depleted sections, centered at x = ± 15 Mm of the initial hux rope, rise in the 
convection zone under buoyancy. The dramatic decrease of pressure in the ambient plasma with 
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altitude gives rise to the vast expansion in the two rising sections of the flux rope as it emerges. 
Fig. [3] illustrates the three-dimensional (3D) structure of the magnetic flux rope in the simulation 
domain after 2.5 hours of evolution with the gray plane showing photospheric vertical magnetic 
field Bz- Each of the two buoyant section rises and expands, bending the initial flux rope and 
forming an M-shaped configuration of magnetic fields extending from the convection zone into 
the corona. At the same time, the mass draining in the buoyant sections down along the field 
lines under gravitational force further enhances the buoyancy in the rising sections. Gravitational 
force of the heavy material keeps the central segment between the two rising sections down in the 
convection zone, while the two buoyant sections emerge. The downward push from the drained 
mass at the central segment between the rising bipoles maintains the configuration of multiplicity 
of bipoles in the domain, and prevents the multiple bipoles from merging into one pair of bipoles. 
Note that in the real solar convection zone, the downward convective flow may make an additional 
contribution to the downward push and facilitate the formation of multiple bipoles. The emerged 
sections expand dramatically due to the stratification of density in the surrounding plasma and 
especially when it approaches the photosphere, density drops by 6 orders of magnitude over a 
thin layer of 2 Mm, driving an enormous expansion of the flux rope in the horizontal direction. 
This expansion gives rise to an unstable situation where high-density material is supported by 
low-density plasma at the top of the flux rope. Under this condition, Rayleigh-Taylor instability 
takes place at the top boundary of the emerging fields and produces the undular structure on the 
field lines as represented by the red lines in Fig. [3j The grey plane of photospheric also shows 
strands of opposite polarities in the c enter of each bipole, man ifesting the effect of Rayleigh-Taylor 
instability on the photospheric plane. Ilsobe et al.l (j2005l . l2006l l have shown that magnetic Rayleigh- 
Taylor instability in simulations of flux emergence into preexisting coronal held also takes place 
in the top-heavy conhguration caused by the expansion of the emerging dome and the draining 
of plasma along the held lines, and produces dense hlaments mimicing Ha arch hlaments, which 
may cause intermittent heating and the hot and cold loops observed in extreme ultraviolet images 
(jYoshimura fc Kurokawalll999l ). 


The initial hux rope has a right-hand twist with q = 0.5, which causes the bipoles to initially 
emerge with its polarity inversion line at an angle to the axis of the hux rope. As the expansion 
of the buoyant emerging bipoles continues, the polarities in each bipole start to shift toward the 
central axis of the initial hux rope. At the center of the domain, the two polarities next to the dip 
in the M-shaped structure are pushed together by the expansion in each of the emerging bipoles, 
as well as the shifting motion toward the central axis, as shown in Fig. 01 These motions drive the 
two opposite polarities together to naturally form a sharp PIT between the the two polarities PI 
and N2 (shown by arrows on Fig. [3]) from the two bipoles. Note that underneath the photosphere, 
the two polarities PI and N2 at the PIT are connected by a subsurface U-loop structure, while 
for each of the bipoles, i.e. Nl-Pl and N2-P2, the opposite polarities are connected by D-loops 
above the surface. At time t = 3.5 hrs, the magnetic fields, shown by the red lines in Fig. 01 expand 
further into the corona, forming overlying coronal loops connecting each of the bipoles. Importantly, 
along the PIT in the center, the magnetic fields from the two opposite polarities PI and N2 are 
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anti-parallel to each other, while their footpoints are pushed together during the emergence and 
expansion. In such a configuration of anti-parallel magnetic helds, strong current sheet builds up 
at the PIL where the two polarities PI and N2 are moving toward each other, as shown by the 
magenta surfaces in Fig. [5] between the red lines connecting Nl-Pl and N2-P2. With the continuing 
convergence of the two polarities during the expansion of the two emerging bipoles, current keeps 
thinning until reconnection takes places at the PIL, between the two nearby polarities PI and N2. 
The three groups of lines in Fig. [6] illustrate the complex coronal magnetic fields in our simulation 
domain. The reconnection forms the large-scale overlying loops connecting the two polarities at 
the two ends of the flux rope N1 and P2, as shown by the magenta lines in Fig. El as well as a 
set of small loops connecting the nearby polarities PI and N2 in the center shown by the blue 
lines. The connectivity of the initial magnetic fields is colored by red, which shows an M-shaped 
structure. By comparing Fig. 0] and El we notice the change of the magnetic connectivity under the 
continuous expansion and emergence. We note that it is the reconnection at the PIL that produces 
the complexity in the structure of magnetic held. At this time, d-spot forms in the center with 
complex magnetic connectivity among the polarities and a compact magnetic structure at the PIL 
(see Fig. E])- 


The evolution of total unsigned hux on the photosphere is shown in Fig. [71 We observe a fast 
growth in the total flux until about time t = 2.6 hrs, when the two Il-shaped structures come into 
contact and collide. Until this time, the expansion and emergence of each structure are independent 
of each other and the magnetic flux simply grows into regions of less magnetic pressure, hlling the 
space between the two buoyant segments rapidly. After time t = 2.6 hrs, the majority of the axial 
hux reaches the photosphere, as shown in Fig. El forming a quadupolar structure on the surface. 
After that, the growth of the total unsigned hux is slowed down, indicated by the less steep gradient 
from time t = 2.6 to 4.5 hrs, when the two polarities PI and N2 collide into each other. During the 
collision, strong current builds up at the sharp PIL, and reconnection takes place between the two 
polarities PI and N2, producing hux cancellation at the PIL. We evaluate the potential magnetic 
energy from the photospheric magnetogram and calculate total free magnetic energy in the corona 
by subtracting the potential energy from the magnetic energy in the coronal helds, as shown by 
the red dashed line in Fig. [71 The free energy keeps increasing with the hux emergence until 3.6 
hrs, when the reconnection at the PIL starts to release the stored free magnetic energy by forming 
shorter loops shown by the blue lines in Fig. El with a total drop of 2.8 x 10^^ ergs in free energy. The 
release of f ree en ergy is comparable with the estimated value of 3.4x10^^ ergs in AR 11158 as in 


Sun et al.l (j2012l ) although the energy release is more gradual in the simulation with no eruptions. 


4. Photospheric Flows and Energy Flux 

At the photosphere, the two buoyant sections of magnetic hux rope centered at x = ± 15 Mm 
hrst emerge as two separate bipoles at time t = 2.5 hrs, as shown in Panel (a) of Fig. El The distinct 
emergence of the two sections can also be viewed in the 3D magnetic structure shown at the same 
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time in Fig. [3l Until then, the emergence of the two fl-shaped structures remains independant 
of each other, with only their footpoints being connected in the convection zone. The blue and 
yellow arrows represent the horizontal velocity in the positive and negative polarities, respectively. 
It is clear that at this time the velocity fields in Nl-Pl and N2-P2 are mostly pointing outward 
from the center of each emerging section, representing the expansion of the emerging domes in a 
stratified atmosphere with dramatic drop in the pressure of the ambient plasma. Due to the initial 
seperation between the two density-depleted sections and the mass draining into the central dip, 
the two polarities PI and N2 in the center remain separated at this time t = 2.5 hrs and unaffected 
by each other with a distance from each other as shown in Fig. O The continuous expansion during 
the emergence in the stratified atmosphere drives the outward motion of each polarities from the 
center of each bipole, as shown by the arrows representing horizontal velocity. Inevitably, the two 
D-shaped emerging domes come into contact with each other and the two opposite polarities PI 
and N2 are pushed together during the expansion, as shown by the velocity arrows in Panel (b) of 
Fig. [8l At the same time, a sharp PIT with high field gradient forms between PI and N2 as they are 
pushed and squeezed into each other, which is found to be the preferential regions of flares and often 
used as an observati o nal predictor of solar erup tive activities (jLeka fc Barnesll2nr),'ll : ISchrijveiiDnnTl : 


Falconer et al.l 120081 : iMason fc Hoeksemal 120111 ). More importantly, as the two polarities cannot 


annihilate with each other during the build up of sharp magnetic gradient, it is necessary that the 
two spots PI and N2 take another curvilinear path around each other during the expansion of the 
two bi poles, whic h is also found in the observations of d-spots formed by collision of non-paired 
spots (|Tand Il983l l. The horizontal velocity in the positive (PI) and negative (N2) polarities, as 
shown by the blue and yellow arrows, respectively, runs almost anti-parallel to each other along the 
PIT, exhibiting a strong shearing patte rn between the two polarities. The velocity shear, driven by 
the spot motion (IZirin fc Liggettlll987l ) ~ flux emergence in our case, persists during the interaction 
of the two polarities PI and N2, as shown by the anti-parallel pattern of velocity arrows in Panel 
(b)-(d) at the PIT. It streches the two sets of magnetic field lines connecting Nl-Pl and N2-P2 
respectively at the PIT in opposite directions and forms a set of highly sheared field lines at the 
interface between PI and N2, as shown in Fig. 01 As discussed in Section[3l reconnection takes place 
at the PIT and for ms a set of complex coronal s tructu res as shown in Fig. [6l The reconnection here, 
similar to that in Ivan Ballegooiien fc Mackavl (j2007l ) where two neigbouring D-loops interact and 
reconnect at the PIT, releases the free magnetic energy by reconfiguring the stretched fields at the 
PIT into more relaxed ones. The continuous shearing motion represented by the arrows adjacent to 
the PIT both in Panel (d) of Fig. [8] and the movie available online then moves the footpoints of the 
post-reconnection low-lying field lines above the PIL in the opposite directions and stretch them 
into elongated and sheared magnetic structures. The highly sheared magnetic conf iguration at the 
PIL is believed to be one of the preceding signatures of great flares in (5-spots (IZirin fc Liggett 


19871 ). during which the free energy is released by relaxing and shortening the magnetic field lines. 


Tana (jl982l l also pointed out that the the polarity axis in (5-spots often deviates from Hale’s law, 
based on sunspot observations from Mount Wilson. In our simulation, as shown by the evolution of 
polarities at the photopshere in Panel (a)-(d) of Fig. [HI we naturally get the inverted polarity of the 
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5-spot PI and N2 with the current configuration of the subsurface flux rope from the collision of 
the two emerging bipoles, as the following polarity of the 5-spot is the leading one in the following 
bipole. During its evolution, however, the rotation motion of the 5-spot drives the two polarities 
PI and N2 in a direction towards obeying the Hale’slaw. The counterclockwise rotation of the 
5-spot is one of the manifestations of the curvilinear paths of the two polarities PI and N2, which 
is ultimately driven by the expansion in each bipole. 


As a result of the cuvilinear path of the spots, the magnetic fields between the two spots take 
the form of highly compressed and streched field lines at the PIL. The blue lines in Fig. Oshow a 
group of sheared reconnected fields among the more relaxed field lines connecting all four polarities 
in the domain shown by the red lines. Fig. [9] shows the strength of the horizontal magnetic field, 
i.e. Bh = + in the 5-spot region outlined by black rectangle in Fig. [H Approaching 

the PIL, the horizontal field strength quickly increases to above 2 kG, as compared to the nearly 
vertical held in the center of each spot shown by the dark blue color. These strong transverse 
magnetic helds with strength comparable to that in the umbrae are conh ned to long, ex t remel y 
narrow channels in 5-spots in Fig. [9] and dehned as “magnetic channels” bv IZirin fc Wand (j 199,11 ). 
In addition, the arrows in Fig. [9] representing the horizontal magnetic helds align almost parallel 
to the PIL, in dicating th at the transverse helds of kG-strength are along the PIL and thus highly 
non-potential (jLowlll982l i. with the free magnetic energy stored which could possibly be released 
in eruptions. 


Besides the shearing motion that forms highly elongated and sheared held lines adjacent to the 
PIL, a converging motion occurs in the outer periphery of the encountering spots, and keeps pushing 
the two spots PI and N2 together. This converging motion results from the expansion and squeeze 
from the emerging bipoles of the writhed emerging rope, as shown by the 3D magnetic structure in 
Fig. m In Fig. [TOl we calculated the Lorentz force in the 5-spot region outlined by the rectangle in 
Panel (d) of Fig. [8l At the outer periphery of the 5-spots PI and N2 away from the PIL, we hnd 
that the horizontal Lorentz forces are pushing the two spots together with a converging pattern, 
consistent with the converging motion as shown in Panel (d) of Fig. [H Meanwhile, at the PIL in 
the middle, apparently there is a strong repelling Lorentz force which pushes the two polarities in a 
direction away from each other. This repelling force naturally develops in such a system when two 
opposite polarities are pushed together with strong current sheet in between, while the converging 
motion in the 5-spots is a manifestation of the expansion of each bipole. The combination of 
repelling force at the PIL and converging motion at the periphery gives the 5-spots the observed 
curvilinear motion around one another. We here point out that this curvilinear motion is due to 
the emergence of the subsurface structure, which pushes a nd closely locks the two polar i ties P I and 
N2 together during the expansion of the two bipoles. In Ivan Ballegooiien &: MartensI (|l989l L the 
presence of converging and shearing motions produces flux cancellation at the neutral line and the 
formation of a flux rope, which is not observed in our simulation here. The absence of such coherent 
a hel ical coronal flux rope may be the result of numerical diffusion (jvan Ballegooiien fc Mackav 


20071 ). given the large cell size of 1237.4 km in the coronal region (see Section [2]). 















To quantify the energy transfer associated with the horizontal and vertical motions of the 
magnetic fields, we calculate the Poynting flux associated with the two velocity components by 
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Fig. 111! shows a zoom-in view of the distribution of the energy fluxes associated with horizontal (a) 
and vertical (b) motions, in the central polarities (same region as Fig. [Toll at time t = 4:40:00. The 
arrows in Panel (a) shows the horizontal flows in the (5-spot and clearly represent a rotation pattern 
in each of the polarities as well as a strong shearing motion adjacent to the PIT shown by the white 
line. Significant energy flux into the corona occurs at the PIT, identified by the two red bands in 
Panel (a) which coincides with the shearing motion. In addition, within each of the polarity, the 
converging motion toward the PIT shown by the arrows also produces a strong energy flux into 
the corona. In Panel (b), the Poynting flux by the vertical motion of the fields also makes a great 
contribution to the energy transport into the corona, with brightennings along the PIT, indicating 
strong energy flux by the emerging of the magnetic fields. The energy flux concentrated at the PIT 
between the two colliding spots transports a great amount of free energy into magnetic fields above 
the PIT. As shown in Fig. [6l the magnetic fields are highly compressed and non-potential above 
the PIT, with a free energy of 8.3 x 10^^ ergs, about 46% of the total magnetic energy of 1.8 x 10^^ 
ergs. This stored free energy may be th e energy source of the great flares that are often observed 
in (5-spot regions (IZirin fc Liggettlll987l ). 


5. Summary and Conclusions 


In an active-region-scale domain, we carry out a simulation of the evolution of a twisted 
magnetic flux rope with two buoyant sections and study the formation of a compact (5-spot like 
structure during the interaction of the two emerging sections as well as the maintenance of such 
compact structure with a realistic size scale. At the beginning, the two buoyant sections rise into the 
atmosphere independently of each other, expanding dramatically in a highly stratified background 
plasma. On the photosphere, the emergence of each section forms one pair of bipoles, each separated 
from the other. The separation in distance decreases as the two sections emerge and expand. The 
fact that the two buoyant sections originate from the same flux rope prev ents the two bipoles from 
escaping from each other by subsurface connection ljToriumi et al.ll2ni4l ). Inevitably the leading 
polarity of the following bipole collides with the following one of the leading bipole, forming a sharp 
PIT with a high magnetic gradient. The continuous expansion during the emergence keeps pushing 
together the two polarities in the center and locks the two spots closely, producing a compact (5-spot 
active region on the photosphere. The two colliding polarities take a curvilinear path around one 
another, producing the observed shearing and rotation motion. At the PIT, the horizontal velocity 
in the two polarities runs anti-parallel with each other. The magnetic fields is highly compressed 
and almost parallel to PIT in the middle of the (5-spot. The sheared transverse magnetic field of kG 
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strength possesses a significant amount of free magnetic energy up to 9.3 x 10^^ ergs, as a potential 
source for eruptions and flares that are often associated with <5-spot active regions. The 5-spot in 
such an act ive-region-scale simulation produces detailed dynamics comparable with observations 
(|Tandll983l i. such as the inverted polarity against Hale’s law, sharp PIL with sheared kG transverse 
magnetic helds, the curvilinear motion of the polarities, as well as the complex coronal magnetic 
connectivity. 


As collision of two bipolar groups is the most popular type of formation of 5-spots (jZirin &: Liggett 


19871 1 , it is crucial to understand mechanisms that closly lock the two bipoles together during their 
emergence and expansion. As has been pointed out bv iToriumi et al.l (|2014l l. the subsurface con¬ 
figuration plays an important role in the formation of a compact photospheric magnetic structure. 
The two emerging bipoles in our simulation, deeply rooted and connected together in the convection 
zone, are not able to separate from each other due to the internal linkage. The emergence and ex¬ 
pansion drive the two spots into a curvilinear path around one another which manifests as rotation 
of the spots and shearing motion at the PIL on the photosphere. Nevertheless, it is the repelling 
force from each 5-spot’s bipolar counterpart that pushes the two polarities in the 5-spot together 
and maintains the integrity of the quadupolar structure during the evolution. In the realistic solar 
convection zone, the subsurface conhguration of multiple buoyantly emerging sections on the same 
flux rope may occur when the fl ux rope encounter s multiple cells of upflowing plasma, which is 


possible in the turb ulent interior (INelson et al.l 


20111 : 


Fan &: Fanejl2014 L Realistic simulations with 


convective motions (jRempel fc Schlichenmaiei:il201ll : iFang et al.ll2012bl ) have to be implemented to 
study the detailed dynamics in the emergence from the interior to the photosphere. In addition, 
in simulations with the effect of Coriolis force on the flux emergence, a retrograde flow naturally 
develops along the flux rope and produces assymetry in the two polarities of each emerging loop, 
with the leading one being more dominant (jFan et al.lll993l : lR,empel fc Cheundl2ni4l ). Therefore, in 
the process where a 5-spot is formed by collison of two bipoles, the fact that the following polarity 
in the 5-spot is the leading polarity of the followin g bipole cou ld also explain the dominance of the 
following polarity oftern observed in such 5-spots (jTandll983l i. 
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Fig. 1.— Initial vertical stratification of pressure (black) and temperature (purple) in the domain. 
The green line shows the distribution of magnetic pressure in a cross-section cut of the flux rope. 



Fig. 2.— Distribution of the initial density along the axis of the flux rope with two buoyant sections 
centered at x = ±15 Mm. 
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Fig. 3.— 3D structure of the M-shaped magnetic field lines (red lines) in the domain at time t = 
2.5 hrs. The plane shows the photospheric magnetogram. 



Fig. 4.— Same as Fig. [3] but at time t = 3.5 hrs. 
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Fig. 5.— The J-spot region at time t = 4.5 hrs with magenta isosurfaces of | J|/|i?| = 100 in regions 
of B/j > 300 G and colored lines representing the magnetic field lines. 



Fig. 6.— Same as Fig. [3] but at time t = 4.5 hrs. The magenta and blue lines are showing the 
overlying and low-lying magnetic field lines after the reconnection, respectively. The red lines show 
the field lines connecting the four polarities. 
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Fig. 7.— Evolution of total unsigned magnetic flux at the photosphere in the computational domain 
and the coronal free magnetic energy. 
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Fig. 8.— Magnetogram at time t = 2.5 (a), 3.0 (b), 3.5 (c) and 4.5 (d) hrs. Blue and yellow arrows 
show the horizontal velocity in positive and negative polarities with \Bz\ > 800 G, respectively. A 
movie showing the evolution of magnetogram with horizontal velocities is available online. 
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Fig. 9.— The horizontal magnetic held with color showing strength and arrows showing directions 
in the area outlined by the rectangle in Panel (d) of Fig. [8] at time t = 4.5 hrs. The black lines 
represent the PILs. 
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Fig. 10.— Vertical magnetic held in the area outlined by the rectangle in Panel (d) of Fig. [8] at 
time t = 4.5 hrs with blue and yellow arrows showing the horizontal Lorentz force in positive and 
negative polarities, respectively. 
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Fig. 11.— Poynting flux associated with the horizontal (a) and vertical (b) motions of the magnetic 
field lines at time t = 4:40:00. The black arrows represent the horizontal flows in regions with \Bz\ > 
800 G and white lines show the PIL. 




